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Abstract: We propose a spectral curve describing torus knots and links in the B-model. In 
particular, the application of the topological recursion to this curve generates all their colored 
HOMFLY invariants. The curve is obtained by exploiting the full S1(2,Z) symmetry of the 
spectral curve of the resolved conifold, and should be regarded as the mirror of the topological 
D-brane associated to torus knots in the large Gopakumar-Vafa duality. Moreover, we derive 
the curve as the large N limit of the matrix model computing torus knot invariants. 
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1. Introduction 

One of the most surprising consequences of the Gopakumar-Vafa duality [21] is that Chern- 
Simons invariants of knots and links in the three-sphere can be described by A-model open 
topological strings on the resolved conifold [44] (see [40] for a recent review). The boundary 
conditions for the open strings are set by a Lagrangian submanifold associated to the knot or 
link. By mirror symmetry, an equivalent description should exist in terms of open strings in the 
B-model, where the boundary conditions are set by holomorphic submanifolds. 

This conjectural equivalence between knot theory and Gromov-Witten theory has been im- 
plemented and tested in detail for the (framed) unknot and the Hopf link. For the framed unknot 
there is a candidate Lagrangian submanifold in the A-model [44]. Open Gromov-Witten invari- 
ants for this geometry can be defined and calculated explicitly by using for example the theory of 
the topological vertex [2], and they agree with the corresponding Chern-Simons invariants (see 
for example [53] for a recent study and references to earlier work). The framed unknot can be 
also studied in the B-model [3, 41]. As usual in local mirror symmetry, the mirror is an algebraic 
curve in C* x C*, and the invariants of the framed unknot can be computed as open topological 
string amplitudes in this geometry using the formalism of [39, 8]. The Hopf link can be also 
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understood in the framework of topological strings and Gromov-Witten theory (see for example 
[24]). 

In spite of all these results, there has been little progress in extending the conjectural equiva- 
lence between knot theory and string theory to other knots and links. There have been important 
indirect tests based on integrality properties (see [40] for a review), but no concrete string theory 
calculation of Chern-Simons invariants of knots and links has been proposed beyond the unknot 
and the Hopf link, even for the trefoil (which is the simplest non-trivial knot). 

In this paper wc make a step to remedy this situation, and we provide a computable, B-model 
description of all torus knots and links. Torus knots and links are very special and simple, but 
they are an important testing ground in knot theory and Chern-Simons theory. As we will see, 
our B-model description does not involve radically new ingredients, but it definitely extends the 
string/knot dictionary beyond the simple examples known so far. 

Our proposal is a simple and natural generalization of [3]. It is known that for B-model 
geometries that describe mirrors of local Calabi-Yau threefolds, and are thus described by a 
mirror Riemann surface, there is an S1(2,Z) action that rotates the B-model open string moduli 
with the reduction of the holomorphic three-form on the spectral curve; this action is a symmetry 
of the closed string sector. For open strings, it was proposed in [3] that the unknot with / units 
of framing is obtained by acting with the S1(2,Z) transformation T-^ on the spectral curve of 
the resolved conifold (here T denotes the standard generator of the modular group), but no 
interpretation was given for a more general modular transformation. As we will show in this 
paper, the B-model geometry corresponding to a (Q, P) torus knot is simply given by a general 
S1(2,Z) transformation of the spectral curve describing the resolved conifold. This proposal 
clarifies the meaning of general symplectic transformations of spectral curves, which play a crucial 
role in the formalism of [19]. Moreover, it is in perfect agreement with the Chern-Simons 
realization of the Verlinde algebra. In this realization, one shows [31] that torus knots are 
related to the (framed) unknot by a general symplectic transformation. Our result can be simply 
stated by saying that the natural Sl(2, Z) action on torus knots in the canonical quantization of 
Chern-Simons theory is equivalent to the Sl(2, Z) reparametrization of the spectral curve. 

In practical terms, the above procedure associates a spectral curve to each torus knot or link. 
Their colored U {N) invariants can then be computed systematically by applying the topological 
recursion of [19] to the spectral curve, exactly as in [8]. In this description, the (P, Q) torus knot 
comes naturally equipped with a fixed framing of QP units, just as in Chern-Simons theory [31]. 
As a spinoff of this study, we obtain a formula for the HOMFLY polynomial of a (Q, P) torus 
knot in terms of g-hypcrgcomctric polynomials and recover the results of [22]. 

Our result for the torus knot spectral curve is very natural, but on top of that we can 
actually derive it. This is because the colored U{N) invariants of torus knots admit a matrix 
integral representation, as first pointed out in the SU{2) case in [34]. The calculation of [34] 
was generalized to U{N) in the unpublished work [38] (sec also [16]), and the matrix integral 
representation was rederived recently in [6, 28] by a direct localization of the path integral. We 
show that the spectral curve of this matrix model agrees with our natural proposal for the B- 
model geometry. Since this curve is a symplectic transformation of the resolved conifold geometry, 
and since symplectic transformations do not change the 1 /N expansion of the partition function 
[19, 20], our result explains the empirical observation of [16, 6] that the partition functions of the 
matrix models for different torus knots are all equal to the partition function of Chern-Simons 
theory on (up to an unimportant framing factor). 

This paper is organized as follows. In Section 2 we review the construction of knot operators 
in Chern-Simons theory, following mainly [31]. In section 3 we focus on the B-model point of 
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view on knot invariants. We briefly review the results of [3] on framed knots, and we show 
that a general S1(2,Z) transformation of the spectral curve provides the needed framework to 
incorporate torus knots. This leads to a spectral curve for torus knots and links, which we 
analyze in detail. We compute some of the invariants with the topological recursion and we show 
that they agree with the knot theory counterpart. Finally, in Section 4 we study the matrix 
model representation of torus knots and we show that it leads to the spectral curve proposed in 
Section 3. We conclude in Section 5 with some implications of our work and prospects for future 
investigations. In the Appendix we derive the loop equations satisfies by the torus knot matrix 
model. 



2. Torus knots in Chern— Simons theory 

First of all, let us fix some notations that will be used in the paper. We will denote by 

Uk = P e^p(f A (2.1) 

the holonomy of the Chern-Simons connection A around an oriented knot /C, and by 

= Tr rUk (2.2) 

the corresponding Wilson loop operator in representation R. Its normalized vev will be denoted 

by 



Wr{IC) = {Ttr \P exp j^^ Aj J . (2.3) 

In the U {N) Chern-Simons theory at level k, these vevs can be calculated in terms of the variables 
[51] 

When i? = □ is the fundamental representation, (2.3) is related to the HOMFLY polynomial 
'^(/C) of the knot /C as [51] 

^□(^) = ,i/2lV-V2 ^(^)- (2-5) 

Finally, we recall as well that the HOMFLY polynomial of a knot 1-L{1C) has the following structure 
(see for example [35]), 

UiK) = Y,Pi{c^)z^\ 2 = _ ^-i/2_ (2.6) 

i>0 

Torus knots and links have a very explicit description [31] in the context of Chern-Simons 
gauge theory [51]. This description makes manifest the natural S1(2,Z) action on the space of 
torus knot operators, and it implements it in the quantum theory. It shows in particular that all 
torus knots can be obtained from the trivial knot or unknot by an Sl(2, Z) transformation. We 
now review the construction of torus knot operators in Chern-Simons theory, referring to [31] 
for more details. 

Chern-Simons theory with level k and gauge group SU {N) can be canonically quantized on 
three-manifolds of the form S x M, where S is a Riemann surface [51]. The resulting Hilbert 
spaces can be identified with the space of conformal blocks of the U{N) Wess-Zumino-Witten 
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Figure 1: The path integral over a solid torus with the insertion of the Wilson line in representation R 
gives the wavefunction defined in (2.10). 



theory at level k on S. When S = has genus one, the corresponding wavefunctions can be 
explicitly constructed in terms of theta functions on the torus [9, 34, 18, 5]. The relevant theta 
functions are defined as 



0z,p(-^,a) = 



exp 



ivrr/ i^ + j^ + 2vri/ + j 



(2.7) 



where r is the modular parameter of the torus T^, Ar is the root lattice of SU{N), a G Ar C 
and p G A^, the weight lattice. Out of these theta functions we define the function 



V'i,p(T, a) = exp (j^^^<^^^ ^i,p{^i a)- (2i 



Notice that, under a modular 5'-transformation, a transforms as 

a a/r. (2.9) 

A basis for the Hilbert space of Chern-Simons theory on the torus is given by the Weyl antisym- 
metrization of these functions, 

A«,p('r,a)= ^ e(u;)V'i,^(p)(r,a), (2.10) 

where W is the Weyl group of SU{N), and 

l = k + N. (2.11) 

The only independent wavefunctions obtained in this way are the ones where p is in the fundamen- 
tal chamber Ti, and they are in one-to-one correspondence with the integrable representations of 
the affine Lie algebra associated to SU{N) with level k. We recall that the fundamental chamber 
J^i is given by A„//Ar, modded out by the action of the Weyl group. For example, in SU{N) a 
weight p = X][=iPiAi is in I'l if 

r 

Y^Pi<l, and Pi > 0, z = 1, • • • , r (2-12) 

i=l 

where r = — 1 is the rank of the gauge group. The wavefunctions (2.10), where p G J-'i, span 
the Hilbert space ?^;(T^) associated to Chern-Simons theory on T^. 
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The state described by the wavefunction Xi^p has a very simple representation in terms of 
path integrals in Chern-Simons gauge theory [51]. Let us write 

P = P + Ar, (2.13) 

where p is the Weyl vector and Aji is the highest weight associated to a representation R. Let 
us consider the path integral of Chern-Simons gauge theory on a solid torus M^i with boundary 
dMf2 = T^, and let us insert a circular Wilson line 

y^(i,o) = Tri? I P exp ^ A] (2.14) 



•^1,0 



along the non-contractible cycle ICifi of the solid torus (see Fig. 2). This produces a wavefunction 
'^{A), where A is a gauge field on T^. Let us now denote by uj{z) the normalized holomorphic 
Abelian differential on the torus, and let 



R 



H = Y,HiXi (2.15) 

i=l 

where Aj are the Cartan matrices and fundamental weights of SU{N), respectively. A gauge 
field on the torus can be parametrized as 

Az = {UaU)"^ dz {UaU) , A^ = {UaU)"^ 0^ {UaU) , (2.16) 

where 

u:T^^ SU{Nf (2.17) 
is a single-valued map taking values in the complexification of the gauge group, and 

Ua = exp ( / io(z')a ■ H - / uj(z')a ■ H] . (2.18) 

\Imr J Imr J J 

In this way, the gauge field is written as a complexified gauge transformation of the complex 
constant connection 



Y^^iHi. (2.19) 



i=l 

After integrating out the non-zero modes of the gauge connection [18, 31], one obtains an effective 
quantum mechanics problem where wavefunctions depend only on a, and they are given precisely 
by (2.10). In particular, the empty solid torus corresponds to the trivial representation with 
Aij = 0, and it is described by the "vacuum" state 

h,p- (2.20) 

We will also represent the wavefunctions (2.10) in ket notation, as \R), and the vacuum state 
(2.20) will be denoted by |0). 

Torus knots can be defined as knots that can be drawn on the surface of a torus without 
self-intersections. They are labelled by two coprime integers (Q, -P), which represent the number 
of times the knot wraps around the two cycles of the torus, and we will denote them by K.q,p- 
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Figure 2: The trefoil knot, shown on the left, is the (2, 3) torus knot. The knot shown on the right is 
the (3,8) torus knot (these figures courtesy of Wikipedia). 

Our knots will be oriented, so the signs of Q, P are relevant. We have a number of obvious 
topological equivalences, namely 

^Q,p - ^P,Q, ^Q,P - ^-Q-P 
If we denote by /C* the mirror image of a knot, we have the property 

^*Q,P = ^Q-P- (2.22) 

This means that, in computing knot invariants of torus knots, we can in principle restrict ourselves 
to knots with, say, P > Q > 0. The invariants of the other torus knots can be computed by 
using the symmetry properties (2.21) as well as the mirror property (2.22), together with the 
transformation rule under mirror reflection 

(YtU^*) {q, c) = {TtUk) {q-\ c-'). (2.23) 

All the knots ICij, with / G Z, are isotopic to the trivial knot or unknot. The simplest non-trivial 
knot, the trefoil knot, is the (2, 3) torus knot. It is depicted, together with the more complicated 
(3, 8) torus knot, in Fig. 2. 

Since torus knots can be put on T^, a ((5,P) torus knot in a representation R should lead to 
a state in ^((T^). As shown in [31], these states can be obtained by acting with a knot operator 

'Wr'^^ : Hi{T^) HiiT^), (2.24) 

on the vacuum state (2.20). If we represent the states as wavefunctions of the form (2.10), torus 
knot operators can be explicitly written as [31] 

where Mr is the space of weights associated to the representation R. In the above description 
the integers (Q, P) do not enter in a manifestly symmetric way, since Q labels the number of 
times the knot wraps the non-contractible cycle of the solid torus, and P labels the number of 
times it wraps the contractible cycle. However, knot invariants computed from this operator are 
symmetric in P, Q; this is in fact a feature of many expressions for quantum invariants of torus 
knots, starting from Jones' computation of their HOMFLY polynomials in [27]. From (2.25) one 
finds, 

MS Mil 
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(2.21) 



ITT/X — ^ V 27riyP • /X 



(2.26) 




Figure 3: The knot operator with labels (1,/) creates a Wilson line which winds once around the 
noncontractible cycle of the solid torus, and / times around the contractible cycle. This corresponds to 
an unknot with / "ribbons," i.e. to an unknot with / units of framing. 



The torus knot operators (2.25) have many interesting properties, described in detail in [31]. 
First of all, we have the property 

"^R^^^P = A/,p+A«, (2.27) 

which is an expected property since the knot /Ci^o leads to the Wilson line depicted in Fig. 2. 
Second, they transform among themselves under the action of the modular group of the torus 
S1(2,Z). One finds [31] 

Mwf'^^M-^ = Wf^'-^^^, M G S1(2,Z), (2.28) 

where {Q, P)M is the natural action by right multiplication. Since the torus knot {Q, P) = (1, 0) 
is the trivial knot or unknot, we conclude that a generic torus knot operator can be obtained by 
acting with an Sl(2, Z) transformation on the trivial knot operator. Indeed, 

Mq,pWJJ'°)m-J, = Wg^'^\ (2.29) 
where Mq^p is the Sl(2, Z) transformation 

Mq,p = ^) (2.30) 

and 7, 5 are integers such that 

Q5-P-f = l. (2.31) 

Since P, Q are coprime, this can be always achieved thanks to Bezout's lemma. 

The final property we will need of the operators (2.25) is that they make it possible to 
compute the vacuum expectation values of Wilson lines associated to torus knots in S'^. In fact, 
to construct a torus knot in §^ we can start with an empty solid torus, act with the torus knot 
operator (2.25) to create a torus knot on its surface, and then glue the resulting geometry to 
another empty solid torus after an S transformation. We conclude that 

Wr{1Cq,p)- ^^1^1^^ , (2.32) 

where we have normalized by the partition function of S^. When performing this computation 
we have to remember that Chern-Simons theory produces invariants of framed knots [51], and 
that a change of framing by / units is implemented as 

Wr (/C) ^ e^^'f^^WR (/C) (2.33) 
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where 

2{k + N) ■ ^^-^^^ 

For knots in §^ there is a standard framing, and as noticed aheady in [31], torus knot operators 
naturally implement a framing of QP units, as compared to the standard framing. For example, 
the knot operator 

wg'-^\ fez, (2.35) 

creates a trivial knot but with / units of framing [31, 26], see Fig. 2. As we will see, the 
same natural framing QP appears in the B-model for torus knots and in the matrix model 
representation obtained in [34, 38, 6]. 

The vcv (2.32) can be computed in various ways, but the most efficient one was presented 
in [46] and makes contact with the general formula for these invariants due to Rosso and Jones 
[45]. One first considers the knot operator 

wSf'°^ (2.36) 

which can be regarded as the trace of the Q-th power of the holonomy around ICifi- It should 
then involve the Adams operation on the representation ring, which expresses a character of the 
Q-th power of a holonomy in terms of other characters, 

chR{U'^) = ^clQchv{U). (2.37) 

V 

Indeed, one finds [46] 

V 

If we introduce the diagonal operator [46] 

Up/q\R) = e'^^'P/^^^\R) (2.39) 
we can write an arbitrary torus knot operator as 

WS?'^) = Dp/qWS?'°)d-;q = r^/QwS?'°V-^/«, (2.40) 

where 

TP/Q=(YI^\ (2.41) 



is a "fractional twist," in the terminology of [42]. The above identity can be interpreted by saying 
that the holonomy creating a (Q,P) torus knot is equivalent to the Q-th power of the holonomy 
of a trivial knot, together with a fractional framing P/Q (implemented by the operator Dp/g). 
As we will see, the same description arises in the B-model description of torus knots. Since 

(2.42) 

the quantum dimension of the representation R, we find from (2.40) that the vev (2.32) is given 

by 

Wr {JCq,p) = J2 <Qe^-«/™-dim,(y). (2.43) 

V 
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This is precisely the formula obtained by Rosso and Jones in [45] (see also [36] for a more 
transparent phrasing). As pointed out above, in this formula the torus knot comes with a 
natural framing of QP units. 

The formalism of torus knot operators can be also used to understand torus links. When Q 
and P are not coprime, we have instead a link Cq^p with L = gcd{Q, P) components. From the 
point of view of the above formalism, the operator creating such a link can be obtained [26, 32] 
by considering the product of L torus knot operators with labels {Q/L,P/L), i.e. 

Wg;)^^ = nwg/^'^/^). (2.44) 

i=i 

As explained in [33] , this can be evaluated by using the fact that the torus knot operators provide 
a representation of the fusion rules of the afhne Lie algebra [31], therefore we can write 

wSg:'.',H.=E<r,...,«.wg'^'^"', (2.45) 

Rs 

where the coefficients in this sum are defined by 

Ri®---®Rl = Y, Kl- ,Ri.Rs (2.46) 

Rs 

and can be regarded as generalized Littlewood-Richardson coefficients. The problem of torus 
links reduces in this way to the problem of torus knots. Notice that in this formalism each 
component of the torus link has a natural framing QP/L"^. 



3. The B-model description of torus knots 
3.1 Preliminaries 

Before discussing the B-model picture, we will recall the standard dictionary relating the correla- 
tors obtained in the knot theory side with the generating functions discussed in the B-model (see, 
for example. Appendix A in [8]). In the knot theory side we consider the generating function 

F(V) = log ZiV), ZiV) = J2 Wr{JC) Tvr V (3.1) 

R 

where V is a C/ (oo) matrix, and we sum over all the irreducible representations R (starting with 
the trivial one). It is often convenient to write the free energy F{V) in terms of connected 
amplitudes in the basis labeled by vectors with nonnegative entries k = (A;i , A;2 , • • • ) . In this 
basis, 

^(^) = J2 -^k'^T^k(V) (3.2) 
k 

where 

oo 

Tk(V) = lliTvWp, = n^^'j"'- (3.3) 

The functional (3.1) has a well-defined genus expansion, 

oo oo 

F{V) = Y^Y.9r'^'^i'\zi,---,ZH). (3.4) 

5=0 h=l 
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In this equation, gg is the string couphng constant (3.8), and we have written 

h 

where myj{z) is the monomial symmetric polynomial in the Zi and is the symmetric group of 
h elements. After setting Zi = p^^, the functionals A^i^\zi, • • • , Zh) are given by 

Ah\pi, ■■■ ,Ph) = J dpi • • • dphWg^hiPi, ■ ■ ■ ,Ph), (3.6) 

where the functionals Wg^h are the ones appearing naturally in the B-model through the topo- 
logical recursion. 

3.2 Symplectic transformations in the resolved conifold 

We now briefly review the B-model description of the framed unknot proposed in [3] . 

According to the Gopakumar-Vafa large N duality and its extension to Wilson loops in [44] , 
knot and link invariants are dual to open topological string amplitudes in the resolved conifold 

0(-l)eO(-l) (3.7) 

with boundary conditions set by Lagrangian A-branes. We recall the basic dictionary of [21]: the 
string coupling constant gs is identified with the renormalized Chern-Simons coupling constant, 

while the Kahler parameter of the resolved conifold is identified with the 't Hooft parameter of 
U{N) Chern-Simons theoy, 

t=P^= gsN. (3.9) 
k + N ^ ^ 

The unknot and the Hopf link correspond to toric A-branes of the type introduced in [4, 44] and 
their Chern-Simons invariants can be computed in the dual A-model picture by using localization 
[29] or the topological vertex [2, 24, 53]. 

By mirror symmetry, there should be a B-model version of the Gopakumar-Vafa large 
duality. We recall (see for example [3] and references therein) that the mirror of a toric Calabi- 
Yau manifolds is described by an algebraic curve in C* x C* (also called spectral curve) of the 
form 

If (e", e") = 0. (3.10) 

We will denote 

= e", V = e'" (3.11) 

The mirrors to the toric branes considered in [4] boil down to points in this curve, and the 
disk amplitude for topological strings is obtained from the function v(u) that solves the equation 
(3.10). Different choices of parametrization of this point lead to different types of D-branes, as we 
will discuss in more detail. According to the conjecture of [37, 8], higher open string amplitudes 
for toric branes can be obtained by applying the topological recursion of [19] to the spectral curve 
(3.10). 
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The mirror of the resolved conifold can be described by the spectral curve (see [3, 8]) 



H{U,V) = V-c-^UV + cU -1 = 0, (3.12) 

where 

c = e*/^ (3.13) 

By mirror symmetry, t corresponds to the Kahler parameter of the resolved conifold. Due to 
the identification in (3.9), the variable c appearing in the spectral curve is identified with the 
Chern-Simons variable introduced in (2.4). The mirror brane to the unknot with zero framing, 
JCifi, is described by a point in this curve, parametrized by U, and the generating function of 
disk amplitudes 

-logV{U) = - log (^^^) = E^Tri^£^ J,=oC/", (3.14) 

^ n>0 

can be interpreted as the generating function of planar one-point correlators for the unknot. 

As pointed out in [3], in writing the mirror curve (3.10) there is an ambiguity in the choice 
of variables given by an Sl(2, Z) transformation, 

X = U'^V^, 

(3.15) 

Y = 

where Q, P, 7, S are the entries of the Sl(2, Z) matrix (2.30). However, only modular transforma- 
tions of the form 

Mi,f = (J () ' / e (3.16) 

were considered in [3] . In the case of the mirror of the resolved conifold they were interpreted 

as adding / units of framing to the unknot. It was argued in [3] that only these transformations 
preserve the geometry of the brane at infinity. The resulting curve can be described as follows. 
We first rescale the variables as 

U,X ^c^U,c^X. (3.17) 



The new curve is defined by. 



X = U 



V 



(3.18) 

1-cf-W 

and as proposed in [39, 8], the topological recursion of [19] applied to this curve gives all the 
Chern-Simons invariants of the framed unknot. 

The general symplectic transformation (3.15) plays a crucial role in the formalism of [19, 8], 
where it describes the groTip of symmetries associated to the closed string amplitudes derived from 
the curve (3.10). It is natural to ask what is the meaning of these, more general transformation. 
In the case of the resolved conifold, and in view of the modular action (2.29) on torus knot 
operators, it is natural to conclude that the transformation associated to the matrix Mq^p leads 
to the mirror brane to a torus knot. We will now give some evidence that this is the case. In the 
next section we will derive this statement from the matrix model representation of torus knot 
invariants. 
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3.3 The spectral curve for torus knots 

Let us look in some more detail to the general modular transformation (3.15). We first redefine 
the X, U variables as 

U c^/^U, X c^X. (3.19) 

This generalizes (3.17) and it will be convenient in order to match the knot theory conventions. 
The first equation in (3.15) reads now 

and it defines a multivalued function 

U = U{X) = X^/'^ + ■■■ (3.21) 
Equivalently, we can define a local coordinate ( in the resulting curve as 



Combining (3.21) with the equation for the resolved conifold (3.12) we obtain a function V = 
V{X). After re-expressing U in terms of X in the second equation of (3.15), and using (2.31), 
we find that the dependence of Y on the new coordinate X is of the form 

logF= ^logX+ilogF(X). (3.23) 

The term log V(X) in this equation has an expansion in fractional powers of X of the form n/Q, 
where n G Z. By comparing (3.22) to (3.18), we conclude that the integer powers of X appearing 
in the expansion of log V{X) are the integer powers of C*^ in the curve (3.18), but with fractional 
framing 

/ = P/Q- (3.24) 

This is precisely the description of {Q,P) torus knots appearing in (2.40)! It suggests that the 
integer powers of X in the expansion of logV{X) encode vevs of torus knot operators. Since 
the first term in (3.23) is not analytic at = 0, we can regard logy (up to a factor of Q) as 
the spectral curve describing torus knots in the B-model. Equivalently, if we want a manifestly 
analytic function of ^ at the origin, as is the case in the context of the matrix model describing 
torus knots, we can consider the spectral curve in the {X, V) variables defined by 



X = U'^ 



1 _ c-P/Q-if/ 



■ 

/ (3.25) 

This curve can be also written as 

Hq^p{X, V) = VP{V - i f - c^~^X{V - c^f = 0. (3.26) 
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Notice that, when Q = 1, P = (i.e. for the unknot with zero framing) we recover the standard 
equation (3.12) for the resolved conifold, and for Q = 1, P = / we recover the curve of the framed 
unknot (3.18). In the curve (3.25), X is the right local variable to expand in order to obtain 
the invariants. The topological recursion of [19], applied to the above curve, leads to generating 
functionals which can be expanded in powers of X^/*^ around X = 0. The coefficients of the 
integer powers of X in these expansions give the quantum invariants of the (Q, P) torus knot, in 
the QP framing. 

When Q and P are not coprime, the above curve describes a torus link with L = gcd{Q, P) 
components. Up to a redefinition of the local variable of the curve, the disk invariants have the 
same information of the disk invariants of the (Q/ L, P/L) torus knot. However, as we will see 
in a moment, the L-point functions obtained from the topological recursion compute invariants 
of the torus link. 

3.4 One-holed invariants 
3.4.1 Disk invariants 

The simplest consequence of the above proposal is that the integer powers of X in the expansion 
of — log V{X) give the invariants 

We will now compute in closed form the generating function — \ogV{X). The equation (3.22) 
defines the local coordinate C as a function of U, and it can be easily inverted (by using for 
example Lagrange inversion) to give, 

oo 

U = J2(^nC, (3.28) 

n=l 

where ai = 1 and 

This is essentially the result obtained in [41], eq. (6.6), in the context of framed knots, but with 
a fractional framing PjQ. Prom this expansion it is easy to obtain 

- log V{X) = ^ Wn(c)X"/«, (3.30) 

n>l 

where 

^^(-) = ti-^r' (D "n ' - j) , (3.31) 

which is again essentially the result obtained in eq. (6.7) of [41]. Integer powers of X corresponds 
to n = Qm, m G N, and we conclude that the the planar limit of (3.27) for the {Q, P) torus knot 
with framing QP should be given by 

mQ , , mQ— 1— ^ 



2Fi{mP,-mQ;mP-mQ + l;c^) . (3.32) 



e=o ^ ^ j=-e+i 

(-l)"^Qc"'^P-Q)(mP - 1)\ 



{mP — mQ)\{mQ)\ 
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This can be verified for the very first values of Q, P. For example, we obtain: 

{TtUk,,)^^^ = (3c^ - + 2) , (3.33) 

(TVZ^K;3.5>g=o = (7c' - 15c' + 10c' - 2) , 

which give the correct result for the genus zero knot invariants. In particular, the above expression 
turns out to be symmetric under the exchange of Q and P, although this is not manifestly so. 

The expression (3.32) can be written in various equivalent ways, and it is closely related to 
a useful knot invariant. Indeed, the vev 

(T^Wk;),=o(c) (3.34) 

is, up to an overall factor of c — c~^, the polynomial po{c^) appearing in the expansion (2.6). 
This polynomial plays a distinguished role in knot theory, and this seems to be closely related 
to the fact that it is the leading term in the large A'^ expansion (this was first pointed out in 
[13]). The polynomial po(c') of torus knots appears in the work of Traczyk [50] on periodicity 
properties of knots, but a closed expression as a function of Q, P does not seem to be available 
in the literature. Using the above results, and performing various simple manipulations, we find 
the following expression, valid for Q, P > 0: 

^ [c?) = c(^-i)(Q-i) ^^ + ^^7^^' ,F,{l-P,l-QA-P-Q; c^) . (3.35) 

Here 2-^1 (a, c; x) is the standard Gauss' hypergeometric function. Of course, since the indices 
are negative, the r.h.s. is a polynomial in In writing (3.35), which is manifestly symmetric 
under the exchange of P and Q, we have implemented two small changes w.r.t. (3.32). First 
of all, invariants of knots in §^ are usually presented in the standard framing, while the results 
obtained for the spectral curve correspond to a torus knot with framing QP. In order to restore 
the standard framing we have to multiply the expression (3.32) by c~^^ . Second, our labeling 
of the torus knot docs not agree with the standard conventions in the literature: what wc call 
the (Q, P) torus knot is usually regarded as a {Q, —P) torus knot. This means that we have to 
apply the mirror transformation (2.23) to our invariant, which implies in particular that 

Po*(c')=P^(c-')- (3.36) 

After implementing these changes, one obtains (3.35) from (3.32). Of course, if {Q,P) are not 
both positive or both negative, we can use (3.36) to compute the invariant. 

The spectral curve (3.25) gives, on top of the invariants of torus knots, information about 
other invariants associated to the torus knot, encoded in the coefficients of the fractional powers 
of X. They correspond to fractional powers of the holonomy around the knot. As we will see in 
Section 4 these invariants have a natural interpretation in the matrix model for torus knots. 

3.4.2 All-genus invariants 

Even more remarkably, the close relation of the invariants of the (P, Q) torus knots to the ones 
of the unknot at fractional framing can be further pushed to derive an all-genus completion of 
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(3.35) in terms of g-hypergeometric polynomials. To see this, notice that one-holed invariants at 
winding number m receive contributions from vevs in hook representations 

{T^^ic) = T.XR(^rn)Wn{lC) = ^ (-lyWR^M (3.37) 

R Rm,3 

where km is the conjugacy class of a length m cycle in Sm, and denotes a hook representation 
with s + 1 rows. For the framed unknot, we have that 

(ttU^^^^) = J2 (-l)^«''^'^'^«'"-dim,(i?„,,). (3.38) 
The quantum dimension of the representation Rm.s can be written as 

„m(m-l)/4-sm,/2 / 1 \ * / <7* \ 

^'-«(^--' = llllra-s-ms/ '" n (l - n - I) ("9) 
where for n G N the g-number [n] and the g-factorial [n] ! are defined as 

[n] = - [n]! = [n][n-l]...[l]. (3.40) 
Upon applying the Cauchy binomial formula 



s=0 

we obtain the finite sum 



J2fg<m+l)/2 + (3.41) 



m 



Tr?/'"- \ r'^e+mf-m( -.vm+e ^ [mf + £ 1]! 



for the framed unknot at winding number m with / units of framing, which can be regarded as 

a g-deformed version of the formulae of [41] for the framed disc. Following exactly the same line 
of reasoning as we did for the planar case, the full un-normalizcd HOMFLY polynomial for {Q, 
P) torus knots is obtained from (3.42) upon sending f ^ P/Q, m ^ Q: 

h [Q~mi]i[p-Q+£]\ 

f—l \QrP-Q\P - 111 

= IF-Qm. + (3.43) 

where the q-analogue of Gauss' hypergeometric function is defined by 

, / , ^ (a;9)n(& : q)n „ ... 
2(t>i[a,b,c;q,z) = y- —z, (3.44 



n=0 



and the g-Pochhammer symbol is given as {a,q)n = [a + n— l]!/[a — 1]! . Upon taking the q ^ I 
limit, we recover (3.32). The natural g-extension of (3.35) leads to the following expression for 
the HOMFLY polynomial of a torus knot, 

,i/2_\-V2 ^(^^.^) = n^MT^ .0, (1 - P, 1 - Q, 1 - P - Q; c^) . (3.45) 
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Again, since P > and Q > 0, the series truncates to a degree d = min(P — 1,Q — 1) polynomial 
in c^. It can be also written as 

1 n^^r- N JP-1)(Q-1) 1 >^ [P + Q-k-l]l 



ndCn p) = c(-P-i)(Q-^) ^ V [P + Q-k-l]\ (_-^^k 2k .3 



which is the result obtained in [22] for the HOMFLY polynomial of a torus knot. 

With (3.45) at hand we can straightforwardly extract the higher genus corrections to (3.35). 
Expanding the g-factorials around q = 1 

[n]\ = (-l)"n!(l - g)>-i"("+i) + i (n^ - n) {q - 1) + (J^{n - 2)(n - l)n + 

+ ^(n - 2)(n - l)(3n - l)n^ {q - if + 0{q - if^ , (3.47) 
we obtain for example the closed expression 

^1 48 dc2 +48(l-c2) dc +48^^^ ^ ^ ^'^^ 

c(^-i)(Q-i)(P + Q-l)! / 2c^(P-l)P(Q-l)Q ,x 
- 24P!Q^ 1^ P + Q-1 2Pi(2-P,2-Q,-P-Q + 2,c) 

+ (P^Q - P2 + PQ2 - PQ - Q2 _ 1) (1 - P, 1 - Q; -P - Q + 1; c2) I . 



(3.48) 



We get for instance 



P^'(c) = ^c2(c' + 10), (3.49) 
pf"^(c) = -^c^(2c^-9), (3.50) 
(c) = (2c4 - 32c2 + 49) . (3.51) 

in complete agreement with explicit computations using the Rosso-Jones formula (2.43). 
3.5 Higher invariants from the spectral curve 

Let us now move to the case of higher invariants by applying the Eynard-Orantin recursion 
[19] to the spectral data (3.25) or (3.15). Let Tq^p ~ CP^ be the projectivization of the affine 
curve (3.26). We will take U as an affine co-ordinate on CP^ and we will keep using X,Y 
for the meromorphic extensions X,Y : Tq^p — >■ CP^ of (3.25); we will finally call {qi} the 
quadratic ramification points of the X — t- CP^ covering map. Following [19], we recursively define 
a doubly infinite sequence of meromorphic differentials u}g^h{Ui,--- , Uh)dX (Ui) . . . dX [11^) G 
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A^'^(Sym'*(rQ_p)), g > 0, h > 1 on the /i*'^ symmetric product of Tq^p as 

ooo,2{UuU2) = B{Ui,U2), (3.52) 
Ug^h+iiUo, Ui...,Uh) = Y^ Res KiUo, U) (Wg_^,h+2{U, U,Ui,..., Uu) 

/=OJC-f/ 

In (3.52), B{z,'w)dzdw is the Bergmann kernel of Yq,p, namely, the unique double differential 
with a double pole at z = and holomorphic elsewhere. Since Fp^Q 2± CP^, it reads simply 

On the r.h.s. of (3.53), U is the conjugate point to U near U = Qi under the X projection (i.e. 
X{U) = X{U), U / U), the recursion kernel K{Ui, U2) is defined as 

= - 2X'(^,) log m)- log y(r7.)' ^'-''^ 

with 7 U J = {?7i, . . . , Uh}, I n J = ^, and denotes omission of the terms {h, I) = (0, 0) and 
{9, J)- 

The identification of (3.25) as the spectral curve associated to {P, Q) torus knots in 
entails the identification of the differentials u)g^h{U{Xi), • • • , U {Xh))dXi . . .dX^ with the con- 
nected generating functions Wg^h{U{Xi), ■ ■ ■ ,U {Xh))dXi . . . dXh of (3.6) for all {g, h) 7^ (0, 2); 
in the exceptional case (g, h) = (0, 2), the annulus function is obtained from the Bergmann kernel 
upon subtraction of the double pole in the X co-ordinate 

Wo,2{X,,X2) = B{X^,X2) - ^x^^x^y (3-56) 

With (3.25) and (3.52)-(3.53) at hand, it is straightforward to apply the topological recursion 
to compute higher invariants for torus knots. For the annulus function we obtain from (3.52) 
and (3.56) that 

^o.2(C/i, U2) = ^ - ^xiU^)-XiU2)r ^^-''^ 

The planar part of connected knot invariants (3.2) in the conjugacy class basis where 
^ - ki = 2 for the annulus function, can then be straightforwardly computed as 



= Resu,=ooResu,=ooX{Uirx{U2rWo,2{Ui, U2)dUidU2, (3.58) 

5=0 
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with ki = Sin + Sim- We find explicitly for Q = 2 



W, 



(c) 



(2,0,0,...) 



=0 4 



(c2-l)Pc2^-4 (c\P + l)-P + l) 



{c'^{P + 1)(P + 2) - 2c2p2 + p2 - 3P + 2) 



(1,1,0,0,...) 



=0 9 



:Pc 



.3P-6 



6 (ci2 _ 1) + 4 - 1)^5 + 24 (c^ + 1) (c^ - 1)' P 



5 d4 



+ (55c^ + 82c^ + 55) (c^ - l)^ + 12 (5c^ + 8c^ + 8c^ + 5) (c^ - l)"" P' 



(3.59) 



3 d2 



+ (31c^ + 440*^ + 48c^ + 44c2 + 3l) (c^ - l)^ P 



(3.60) 



and for Q = 3 



W, 



(c) 

(2,0,0,...) 



3=0 24 



Pc^^-^ (^c^^(P + 1)2(P + 2f{P + 3) - 6c^°P(P + lf{P + 2f 
+ 3c^P{P + 1)^(5P(P + 1) - 4)4c^P (5P'' - 7P2 + 2) + 3c^(P - ifP 
X (5(P - 1)P - 4) - Gc^P (p2 - 3P + 2)^ + (P - 3) (p2 - 3P + 2)M , 



(3.61) 

in agreement with the corresponding knot invariants; notice that the case of torus links is also en- 
compassed as soon as gcd(P, Q) > 1, with the Hopf link invariants appearing as the (P, Q) = (2, 2) 
case. 



To compute higher order generating functions we resort to (3.53). The regular branch points 



are 



9± 



c§-i (±^/^^^{c'^ _ 1) p2 + 2 (c2 + 1) PQ + (c2 - 1) Q2 + (c2 _ 1) p + + 

2Q 

(3.62) 

and as will see they are precisely the ramification points that lie on the physical sheet of the 
spectral curve. For the case ^ = 1, ^ = 1 we obtain 

u^i,i{U) = (c2 - 1) (-c^+^) (c2^/3 - (^81c"-f+^ + 81c^U^ - m^(c\P + 3)(P(2P + 3) 

+ 9) - (P - 3)(P(2P - 3) + 9))c^+^ - QU((?{P + 3)(P(2P + 3) + 9) - (P - 3) 
(P(2P - 3) + 9))c-f'+^ + (c^(P + 3)^ - 2c^ (P^ - 54p2 - 162) + (P - 3)^) 

^12 (^3c¥+^ - [/ (c^(P + 3) - P + 3) c^/^ + ?,cU'^Y^ 

(3.63) 



.2P/3 
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and it is immediate to extract genus one, 1-holed knot invariants as 



where in this case ki = Sin- For example 



= Res[/=oo^(t^)"wi,i(C/)dC/, (3.64) 

9=1 



^(1,0,0,...) 
^(1,0,0,...) 



1 /'„2 i\ „P-1 



W, 



(1,0,0,- 



1) Q = 1 

9=1 ^ ^ 

= ^ (c2 - 1) {c\P + 1)(P(P + 2) - 4) + (1 - P)((P - 2)P - 4)) Q = 2 
= TiiC-^"^ c6(P + 1)(P + 2)(2P(P + 3) - 9) - 3c4p(P + 1)(2P(P + 1) Q = 3 

- 5)3c2(P - 1)P(2(P - 1)P - 5) + (2 - P)(P - 1)(2(P - 3)P - 9) 



(3.65) 

which reproduce (3.48) at fixed Q. Similarly, higher winding invariants can be found to reproduce 
the correct knot invariants. 

4. The matrix model for torus knots 

In this section wc study the matrix model representation for quantum, colored invariants of torus 
knots. We first give a derivation of the matrix model which emphasizes the connection to the 
Rosso-Jones formula (2.43), and then we use standard techniques in matrix models to derive the 
spectral curve describing the planar limit of the invariants. 

4.1 A simple derivation of the matrix model 

The colored quantum invariants of torus knots admit a representation in terms of an integral over 
the Cartan algebra of the corresponding gauge group. Such a representation was first proposed 
for SU{2) in [34], and then extended to simply- laced groups in [38] (see also [16]). More recently, 
the matrix integral for torus knots was derived by localization of the Chern-Simons path integral 
[6] (another localization procedure which leads to the same result has been recently proposed in 
[28]). 

The result obtained in these papers reads, for any simply-laced group G, 

Wr{}Cq,p) = ^ I dne-"V2§. -Q ^^^^^ ^ .^y, ^ch^(e"). (4.1) 
Zq,pJ 2P 2Q 

In this equation, 

Zq,p = fdu e-«V29. Yl 4 sinh ^ sinh ^ , (4.2) 
the coupling constant Qs is 

27ri 

9s = PQQs, 9s = 7— — , (4.3) 
k + y 

y is the dual Coxeter number of G, and u is an element in Kyj ® R. a > are the positive roots. 
Notice that, although Q, P are a priori integer numbers, the integral formula above makes sense 
for any Q, P. 
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The easiest way to prove (4.1) is by direct calculation. In order to do that, we first calculate 
the integral 

f dne-"'/2^« n ^si^h ^ si^^ ^chR(e"), gs = fgs, (4.4) 



a>Q 



where / is arbitrary. We will also denote 

Z,j = I dne-V2^. -Q 4,i,h sinh ^ ' " 



a>0 



2/ 



(4.5) 



Let A/j be the highest weight associated to the representation R. Weyl's denominator formula 
and Weyl's formula for the character give, 



Jj2sinh^ = e{w)e 



w{p)-u 



a>0 



wew 



chi^(e") - ^"'^^ 



(4.6) 



and the integral (4.4) becomes a sum of Gaussians, 



V eiw)e{w') f duex.pl + w{p + Ar)-u + w'{p) ■ u/ f\ (4.7) 



Up to an overall factor which is independent of Ar (and which will drop after normalizing by 
Zij), this equals 



exp 



9sf 



e(w) exp {gsp ■ w{p + Ar)) 



(4.8) 



We then obtain 



— / dne-"'/2fl. rr 4sinh — sinh ^chR(e'^ 
Zif J -L-L 2 2/ ^ 

^'J a>0 ■' 



exp 



9sf 



((A. + pf - pA] E..W^H exp(,.p-^P + A,)) ^ ,2./..,.^^) 



(4.9) 



With this result, it is trivial to evaluate (4.1). The change of variables u = Qx leads to 

Wr{}Cq,p) = ^ / dxe--'/25^ n4smli^sinh^ch;^(e«-), (4.10) 



where 

/ = P/Q- (4.11) 

We can now expand cliR^e^^) by using Adams' operation (2.37). The resulting sum can be 
evaluated by using (4.9), and one obtains 

Wr{K:q,p) = Y <Qe2-«/^'^-dim,(F), (4.12) 

V 
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which is exactly (2.43)^. Therefore, (4.1) is manifestly equal to the knot theory result, and in 
particular to the formula of Rosso and Jones for torus knots invariants. Notice that this matrix 
integral representation also comes with the natural framing QP for the {Q,P) torus knot. A 
similar calculation for Zq^p shows that, up to an overall framing factor of the form 

(4.13) 

the partition function (4.2) is independent of Q, P. This can be also deduced from the calculation 
in [16]. 

We also note that there is an obvious generalization of the matrix model representation (4.1) 
to the torus link {Q, P), given by 

>»'(....H.,(£«.p) = J d--'''^'- n/-"i'^»-i'i3^n*%(="-). (4.14) 

Since (4.1) can be calculated exactly at finite N, and the result is identical to (2.43), what 
is the main interest of such a matrix model representation? As in the case of the Chern-Simons 
partition function on S^, it makes possible to extract a geometric, large N limit of the torus knot 
correlation functions, as we will now see. The fact that Zq^p is independent of Q, P up to a 
framing factor strongly suggests that the spectral curves for different Q, P should be symplectic 
transforms of each other. We will verify this and derive in this way the results proposed in section 
3. 



9s ( P .Q \ 2 



4.2 Saddle— point equation 

We will now solve the matrix model (4.1), for the gauge group U{N), and at large N. The first 
step is to derive the saddle-point equations governing the planar limit. An alternative route, 
which provides of course much more information, is to write full loop equations of the matrix 
model and then specialize them to the planar part. This is presented in the Appendix. 



As in [48, 49], we first perform the change of variables 



Ui = PQ log Xi, 



(4.15) 



which leads to 



TV 



i=l 



i<j 



Y[ dme- n ( 2 si^h ^ ) ( 2 sinh 



2P 



Ui — Uj 



2Q 



N 



(4.16) 



i=l i<j 



N 



X JJexp 



i=l 



PQ 



(logXj)" 



P + Q 



(iV- 1) + 1 logx. 



direct calculation of the integral (4.1) is presented in [49] by using the formalism of biorthogonal polynomials. 
The result seems to agree with the above calculations, but the framing factor is not clearly identified. 
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The matrix integral can then be written as 

r ^ 

z = {PQf / n n (-? - (-r - ) 
•'^+ i=i i<j ^ ^ 

N 

X exp 



(4.17) 



i=l 



^(logxO' - (^(iV - 1) + 1 ) logx, 



where the variables Xi = q^^I^^ are thought of as the eigenvalues of a hermitian matrix M of 
size N X N, with only real positive eigenvalues {xi € 



Define now the resolvent 



G{x) = Tt 



N 



X — M ^ — ' X — Xi 
1=1 



(4.18) 



Our observables are expectation values of product of resolvents, and their expansion into powers 
of gs ■ The 1-point function is 



Wix) = {Gix)) = ^gl3-'W,ix) 

5=0 



(4.19) 



and its leading term Wq (x) is called the spectral curve of the matrix model. The 2-point function 
is 

oo 

W2{xi,X2) = {G{xi) G{x2)f^ = {G{xi) G{x2)) - {G{xi)) (G(x2)) = ^ g^s' ^^(a;!, X2) (4.20) 

g=0 

and similarly, the connected n-point correlation function is the cumulant of the expectation value 
of the product of n resolvents 



W„(X1, ...,Xn) = {G{xi) Gix2) . . . GiXn))^'^ = ^ 53'"+"-' W,(xi, . . . , X J 

g=0 



We will denote in the following 



/27ri 

- = exp 



and the 't Hooft parameter of the matrix model is, as usual, 

t = gsN. 



(4.21) 



(4.22) 



(4.23) 



The saddle-point equations for the matrix integral (4.17) are simply 



Pxf 



p p 
xf - xj 



+ 



Qx 



Q 



Q Q 
xY -xJ 



^logx, + ^(iV-l) + l. 
9s ^ 



(4.24) 
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If we use the identity 



= E (4-25) 



— ^ X — bj 

fe=o ^ 



we can write the first term in (4.24) as 



rj-. rr, . ^ ^ ^ rp. l,)~^Q'r ' ^. — fcP/r. . ' (4.26) 

j^j ik=l ^ ^3 k=l j^i ^ ^3 



But 



N 

y ""^Tr^ = y -Tn = -WoixiUJ^'^) (4.27) 

^-^ Xi — UJ "'^Xj Xi — Ul '^^Xj 1 — U '^^ Qs I — U '^^ 

3+i ■' j=l ■' 



and the equation of motion reads, at leading order in l/iV, 
PA-O 

PQlogx + i = Wo(x + 10) + Wo(a; - iO) + ^ W^ixJ"^) + ^ WQixJ"^). (4.28) 

fc=i fe=i 

This is the equation we will now solve. 

4.3 Solving the saddle point equations 

The resolvent VFo(2;) is analytic in C \ C, where C is a finite set of cuts in the complex plane. It 
satisfies 

lim Wo(x) = 0, (4.29) 



and 



lim Wq{x) = t. (4.30) 

a;— >oo 



We now write the exponentiated version of the resolvent as 

y = (7oxe-'^^°(^) (4.31) 



where 

Co = -eW\ (4.32) 
y is analytic in C \ C and satisfies the equation, 

P-i Q-i , 

y{x + iO) n yi^^""^) n = ( _-nV ^"^'^^^ 

k=i k=i 

as well as the boundary conditions 

y{x) ~ Cox, X ^ 0, (4.34) 

and 

y{x) ~ Cq-^x, X — >■ oo. (4.35) 
Notice that y vanishes only at x = and diverges only at x = oo. 
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We now introduce the P + Q functions 



Fkix) = n yi^ io'^'^^^^), 0<k<Q-l, 



=0 



FQ+1^^) = H ^^^^kPMQY 0<1<P-1. 

k=0 ' 



(4.36) 



If we assume that y{x) has a single-cut C on an interval [a, 6], then F]^ has cuts through the 
rotations of this cut by angles which are integer multiples of 21^ /PQ. 

li{)<k<Q — \ and < Z < P - 1, both and Fqj^i have a cut across oj^^^^^C, and 
according to (4.33), under crossing the cut we have 

Ffc(x-iO) = FQ+;(x + iO). (4.37) 

This implies that the function 

P+Q-l 

S{xJ)= n U-Fk{x)) (4.38) 
fc=o 

has no cut at all in the complex plane: 

<S(x + iO,/)=<S(x-iO,/). (4.39) 

Its only singularities may occur when y = oo or when y = (indeed y appears in the denominator 
in (4.36)), and thus the only singularities are poles at x = or x = oo. If we write 

P+Q 

S{x, /) = ^ (-1)^= Sk{x) f+^-^ (4.40) 
fe=o 

then each Sk{x) is a Laurent polynomial of x. Besides, it is clear that 

5(a;x,/) = 5(x,/), (4.41) 

therefore each Sk{x) is in fact a Laurent polynomial in the variable x^*^. We clearly have 

P+Q-i 

5o(x) = l, Sp+q{x)= n Fk{x) = h (4.42) 

fe=o 

as well as the boundary conditions 

o<*<g-i ^tj~-<-;!:c^^°-f; ("3) 

Fk{x) ~ - (-1)-^ OJ'^^ Cq X^ , X ^ DO, 
0<Z<P 1 Pq+/(x)~-(-1)Q^-'«'Co-«x-Q, x^O, 
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This shows that the symmetric functions of the F^s must satisfy 



ri<A:<P-l, Sk{x) = 0{jr'"^), 

x^O lk = P, Sk{x) = (-1)^+Q Cq-^^ x-^Q (1 + 0{x)) 

(p + l<k<P + Q-l, Sk{x) = O {x-P^x^^-P'^P) , 



(4.45) 



and 




P+Q-i 



Six,f) = r^Q + l + {-lfC-''^x-^Qf^ + i-lfC-''^x^^f+ J2 ^fc/' (4-47) 



fe=i 

where are constants. The functions / = Fk{x) and / = Fq+i{x) must all obey this algebraic 
relationship between x^'^ and /: 



We still have to determine the coefficients Sk- 
4.4 Derivation of the spectral curve 

Our matrix model (4.1) can be regarded as a perturbation of a Gaussian matrix integral, and thus 
the resolvent should have only one cut, i.e. the spectral curve must be rational. This determines 
the coefficients s^. 

Saying that an algebraic equation S^x^'^ , /) = is rational means that there exists a rational 
parametrization of the solution. Since the equation is of degree 2 in x^^, this means that, for 
each /, we have two possible values for x^'^, i.e. two points on the spectral curve. In other 
words, / is a rational function of degree 2 of an auxiliary parametric variable which we call V 
(later we shall see that it indeed coincides with the function V{x) defined in (3.25)). Upon a 
Moebius change of variable on V, we can always fix 3 points, and assume that / has a pole at 
V = oo and at F = 1, and a zero at F = 0, i.e. we write it 



where the location of the second zero is to be determined later, but will eventually agree with 
the value given by the definition (3.13). Since the equation is of degree P + Q in /, this means 



where 'R-p+Q is a rational function with P + Q poles. Moreover, the behavior at a; — > (i.e. at 
y ^ or y — > c^) can be of the form 



5(x,/) = 0. 



(4.48) 



f = AV 



l-V 



(4.49) 




f = Fk = Oix^) = 0{x^). 



(4.51) 
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Since / is a rational function, it cannot behave like a fractional power, therefore x^^ must have 
a zero of an order which is a multiple of Q, let us say at V = c^. The behavior of 

/ = Fq+i = 0{x-Q) (4.52) 

implies that x^^ must have a pole of an order which is a multiple of P, let us say at y = 0. 
Similarly, the behaviors at a; — >^ oo, i.e. V ^ 1 or V ^ oo, imply that a;^*^ has a pole of an order 
which is a multiple of Q and a zero of an order which is a multiple of P. Since the total degree 
of x^^ is P + Q, this means that the orders of the poles and zeroes must be exactly P and Q, 
respectively, and there are no other possible poles and zeroes. We have then obtained that 

x^^ = BV-^ (4.53) 

Matching the behaviors of (4.43) and (4.44) gives the values of the coefficients A, c, B: 

A = -l, c = e*/2, B = c^+^, (4.54) 
and in particular identifies c = e*/^ with the variable introduced in (3.13). We finally obtain. 



f = -v 



1 _ c-2y\ Q 



1- V 



(4.55) 



1 - V 

Notice that the relation between X = x~^'^ and V is precisely (3.26). 

To complete the derivation of our spectral curve, let us recall the relationship between / and 
the resolvent. We have, by definition of the resolvent, 

/ 1 oo cx) / N \ 

W{x) = {G{x)) = /Tr-^) =^x-^ (TrM'^)=^x-' /^e'^-^/^^) (4.56) 

^ ' fe=0 fe=0 \i=l / 

and 

P-l P-l , V , p 



1=0 1=0 ^ ' ^ 

oo / JV \ 



(4.57) 



k=0 \ i=l 



In the planar limit we obtain. 



^ ' g=0 ^ ,=0 



P-l 



PQ , y{JQx) 



In 



Q Fo(.x) 



P + Q {-if-^c^xP 

1 (-Fo(x))Q 



P + Q {-Co)PQ xPQ ■ 
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In other words, the resolvent of is (up to trivial terms), the log of F(){x), which is one branch 
of the algebraic function /. From our explicit solution (4.55) we find that 

therefore 

fe=0 \i=l / ^ 

Wc have then proved that t — lnV{x) is the resolvent of . It allows to compute expectation 
values of traces of powers of M which are multiples of P. Our derivation also shows that the X 
in (3.30) is 

X = x-^^. (4.61) 

Since the relation between X = x~^^ and V is precisely (3.26), we have derived the torus knot 
spectral curve from the matrix model. 

We can also deduce from this derivation an interpretation for the coefficients of the fractional 
powers of X appearing in the calculations of section 3: they compute the correlators of the more 
general operators 

chfl (|e""/«) (4.62) 

in the matrix model, which should correspond to "fractional holonomies" around torus knots 
in Chern-Simons theory. Finally, we slioTild mention that the same method used to derive the 
spectral curve makes possible in principe to compute the 2-point function (3.56), and to prove 
the topological recursion (3.52). 



Iny(x). 



(4.60) 



5. Conclusions and prospects for future work 

In this paper we have proposed and derived a spectral curve describing torus knots and links in 
the B-model. The curve turns out to be a natural generalization of [3]: one has just to consider 
the full Sl(2, Z) group acting on the standard curve describing the resolved conifold. Our result 
fits in very nicely with the construction of torus knot operators in Chern-Simons gauge theory, 
and with the matrix model representation of their quantum invariants. 

As we mentioned in the introduction, the ingredients we use to deal with torus knots are 
the same ones that were used to deal with the framed unknot. Going beyond torus knots in 
the context of topological string theory (for example, the figure-eight knot) will probably require 
qualitatively new ingredients in the construction, but this is already the case in Chern-Simons 
gauge theory, where the colored invariants of generic knots involve the quantum 6j coefficients 
[52, 14]. We hope that the results for torus knots obtained in this paper will be useful to address 
the more general case. 

The structure we have found for the invariants of torus knots should have an A-modcl coun- 
terpart. In the A-model, framing arises as an ambiguity associated to the choice of localization 
action in the moduli space of maps with boundaries [29], and the open string invariants depend 
naturally on an integer parametrizing this ambiguity. Our analysis indicates that there should 
be a two-parameter family of open string invariants generalizing the computations made for the 
unknot. These open string invariants can be in principle computed in terms of intersection theory 
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on the moduli space of Riemann surfaces. As an example of this, consider the coefficient of the 
highest power of c (which is P + Q) in the HOMFLY invariant of the (Q, P) torus knot (3.31) 
with framing QP. This coefficient, which we call apj^Q{q), can be expanded in a power series in 
Qs'- 

ap+Q{q) = Y,af^Q9l'-'- (5-1) 
It is easy to see from the results of [29, 41] that the coefficients ctp+g are given by: 

(,) _ (-l)^Q rrr,-| PW.c. f c,iE^i^))c9i^^ii-P/Q - l)n))c,(EV(Pn/g)) 

In this formula, as in [29, 41], Mg^i is the Dclignc-Mumford compactification of the moduli 
space of genus g, 1-pointed Riemann surfaces, ^pl is the Chern class of the tautological line 
bundle Li — >■ Mg^i, E is the Hodge bundle over Mg^i, and we denote 

9 

Cg{E''{u)) = Y,C9-i{^'')u\ (5.3) 

1=0 

Although we have written down an example with h = 1, the generalization to higher h invariants, 
in the spirit of [41], is immediate. Perhaps these formulae can lead to an explicit A-modcl 
description of torus knot invariants, and in particular shed some light on the proposals for 
the corresponding Lagrangian submanifolds made in [33, 47, 30]. Notice that, according to our 
description, the A-model invariants should involve some sort of fractional framing. Such framings, 
in the context of A-model localization, have been considered in [15, 10, 11]. Of course, it should 
be possible to implement the general symplectic transformation we are considering directly in 
the topological vertex. 

Although in this paper we have focused on the spectral curve of the resolved conifold, one can 

consider general Sl(2, Z) transformations of open string amplitudes defined by arbitrary spectral 
curves. In some cases, these transformations have a knot theory interpretation. For example, the 
outer brane in local x describes the unknot in L[2, 1) = MP^ [8], and its general modular 
transformations should describe torus knots in this manifold. It would be also interesting to see 
what is the relation between the approach to torus knots in this paper and the recent work based 
on Hilbert schemes of sing ularities in [43]. 

The topological recursion of [19], which computes open and closed topological string ampli- 
tudes in toric Calabi-Yau manifolds, might have a generalization which gives the mirror of the 
refined topological vertex [25] (recent work in this direction can be found in [1]). If the only data 
entering in this generalization turn out to be the same ones appearing in the original recursion 
(i.e., if the refinement only requires the knowledge of the spectral curve and of the natural differ- 
ential on it, as it happens for example in the P deformation [12]), then one should be able to use 
our spectral curve (3.25) to refine the colored HOMFLY polynomial of torus knots. The resulting 
refinement should provide interesting information on the Khovanov homology of torus knots and 
might lead to a computation of their "superpolynomial" [23, 17], as well as of its generalizations 
to higher representations. 

Finally, the techniques developed in this paper to analyze the matrix model for torus knots 
will probably be very useful in order to understand the large N limit of the more general matrix 
models for Seifert spheres introduced in [37]. Such a large A'' limit would give a way to derive 
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the dual string geometries to Chern-Simons theory in more general rational homology spheres 
-a dual which has remained elusive so far. 
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A. Loop equations 

In this Appendix we derive the loop equations for the matrix integral (4.17). In the following it 
is useful to notice that the resolvent G{x) defined in (4.18) satisfies 



As it is well-known, the method of loop equations consists in observing that an integral is 
unchanged under change of variables. In our case we shall perform the infinitesimal change of 
variable Xi Xi + eSxi -|- O(e^) where 




(A.l) 



Q-lP-l 



,aP+bQ 




(A.2) 



The loop equation, which computes the term of order 1 in e, can be written as 




(A.3) 



where 




(A.4) 



We thus have to compute 





p-1 



(A.5) 




Q-l Q-1 P-1 



a=0 a'=0 b=0 
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We also have 



\ ^ CUXj _ \ ^ \ ^ \ ^ 

^ dxi ~ 2^ 2^ 2^ dxi \ oj"-P+^x - Xi 

i ' 6=0 a=0 i 



= EEE 



6=0 a=0 i ^ 



— X 



2 \ - \ - \ - ^ 

2^ 2^ 2^ (uj'^P+bQx - Xi) 
6=0 a=0 i ^ ^' 

Z-,- Z-,- ^z-, .hO^\0 ^Q\2 2-^ 2-^ 



2 

a=0 6=0 6'=0 

^-^""-^ I ' P + Q-2\ oj'^P+'^Qx 



0=0 6=0 



I.e. 



(A.6) 



^ g (c^ X) _ G{u^''+'Qx). 

6=0 i Vl^ ■^z' -^i .1 a=0 6=0 

The loop equation then gives 

0=0 a'=0 6=0 

+ ^ E E E (gK^+^^x) G(a;«^+^'«x)) (A.7) 

a=0 6=0 6'=0 
Q-iP-1 / . 

EE{Ehn-o + 



^ EEE(^(-"''^'''-)G(-"''''"'''-)) + ^ ^^^(G(a;«^+^«.)GK^+^'«x); 

0=0 o'=0 6=0 a=0 6=0 6'=0 

0=0 6=0 ^ 2 y \ 

Z-,- Z-,- \ 2-^ ,,,aP+bQ^ 



0=0 6=0 ^ \ % 



(A.8) 
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This is then our main loop equation: 



2 ^ ^ ^ \ / 2 

a=0 a'=0 6=0 a=0 6=0 6'=0 

a=0 6=0 ^ ^ 

_^'y^'PQ^aP^,Q^ /y^ ln(a;-^+''Qx)-lnxA 



(A.9) 

We deduce that the spectral curve W^ix) satisfies 

\ If If eVoK^+'''^x)W^oK'^+'''5x) + ^ If Xiif ^o(c^'^^+''^x)W^oK^+''«x) 

a=0 o'=0 6=0 a=0 6=0 6'=0 

= gs If If In (a;'^^+^«x) + ^^^iv) W'oK^+'^a;) - P{x) 

a=0 6=0 ^ 

(A.IO) 

where 

a=0 6=0 \ j * / 

We will now assume that Wq{x) has only one cut C in the complex plane. P{x) has no 
discontinuity through C, 

P(x + iO) = P{x - iO), (A.12) 
and the functions Wo(ci;"^+'''^x) with (a, 6) / (0, 0) have no cut either, so it follows that 

Q-i p-l 
Wo(x + i0)2 + ^ Wo(x + iO)Wo(w"'^x) + ^ Wq{:x + iO)VFo(a;''''3x) 

a'=l 6'=1 



PQ Inx + ) W^o(a; + iO) 



p-l 



(A.13) 



Wo(x - i0)2 + ^ Wo(a; - iO)Wo(w"'^x) + ^ VFo(x + \^)Wq{J''^x) 



a'=\ 6'=1 



- ( PQlna;+ ^^^ t ) Wo(x-iO), 



i.e., if we put all terms in the left hand side and divide by W{x + iO) — Wix — iO), we get the 
equation (4.28) which we derived with the saddle-point method. 
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